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Abstract 

Nitrogen doping in graphene has important implications in graphene-based devices and catalysts. 
We have performed the density functional theory calculations to study the electronic structures of 
N-doped graphene with vacancies and Stone- Wales defect. Our results show that monovacancies 
in graphene act as hole dopants and that two substitutional N dopants are needed to compensate 
for the hole introduced by a monovacancy. On the other hand, divacancy does not produce any 
free carriers. Interestingly, a single N dopant at divacancy acts as an acceptor rather than a 
donor. The interference between native point defect and N dopant strongly modifies the role of 
N doping regarding the free carrier production in the bulk tt bands. For some of the defects and 
N dopant-defect complexes, localized defect n states are partially occupied. Discussion on the 
possibility of spin polarization in such cases is given. We also present qualitative arguments on 
the electronic structures based on the local bond picture. We have analyzed the ls-related x-ray 
photoemission and adsorption spectroscopy spectra of N dopants at vacancies and Stone- Wales 
defect in connection with the experimental ones. We also discuss characteristic scanning tunneling 
microscope (STM) images originating from the electronic and structural modifications by the N 
dopant-defect complexes. STM imaging for small negative bias voltage will provide important 
information about possible active sites for oxygen reduction reaction. 

PACS numbers: 68.43.Bc, 31.15.A-, 82.45. Jn, 82.65.+r 
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I. INTRODUCTION 



Graphene, because of its unique electronic structure, attracts strong attention in the ba- 
sic science field as well as in the application-oriented field. The direct way of controlling 
its physical and chemical properties is to dope different elements. For example, nitrogen- 
doped graphene (N-graphene) has shown wide applications in n-type graphene-based field- 
effect transistors,'^' electrochemical biosensors,^ anodes of lithium-ion batteries,^ catal- 
ysis for various chemical reactions,^ and so on. Among them, we pay particular attention 
to the carbon alloy catalysts (CACs) which are regarded as among the strong candidates 
of Pt-substitute electro-catalysts for oxygen reduction reaction (ORR) at the cathode of 
a fuel cellP^H^ The basic structural components of CACs are multilayered nanographene 
(nanographite) including carbon nanotubes. Although intensive study has been made, the 
reaction site and reaction process are not yet well understood. The main motivation of the 
present study is to clarify the basic aspects related to N doping to graphene so that we can 
obtain some insights into the possible reaction sites of ORR. However, as the basic structural 
and electronic properties of N-graphene form a general background of physics and chemistry 
of graphene, we present the results of our extensive study of electronic structures for various 
local structures of N-graphene in a general form so that researchers working on graphene 
with various aims may utilize the information given below. 

Recent experimental studies^- 4 ' 18 1 20 1 suggest that some mutual effects between structural 
defects and N doping are present during the incorporation of N into graphene. In our 
previous paper,^ from the total energy calculations based on density functional theory 
(DFT) we found that the native point defect (NPD) and N dopant attract each other, 
i.e., they display cooperative effect, suggesting that N dopants can prompt the creation 
of point defects and vice versa. Meanwhile, we have analyzed the energetic stability of 
substitutional N dopants in graphene with vacancies and Stone- Wales (SW) defect. In the 
controllable N doping of graphene by NH3 annealing after ion irradiatiorP or by NH3 plasma 
exposure™ stepwise carrier doping from p-type to n-type has been observed in the measured 
conductivity of N-graphene. 

In the present paper, we concentrate our analysis on how the electronic structure of 
graphene is modified by NPDs and N dopants. Defect states appear at each of monovacancy 
(MV), divacancy (DV) and SW defect in the undoped graphene. The characteristic features 
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of defect states are analyzed in detail for each type of defect. Then the effect of N doping 
on the electronic structures is discussed. We have found several interesting consequences of 
the coexistence of NPDs and N dopants concerning the carrier doping in the bulk 7r bands. 
Although N dopant is generally regarded as electron donor, the N-NPD interaction greatly 
modifies such a naive picture. Several examples of non-trivial behaviors of N dopants are 
presented below. For some of NPDs and N dopant-NPD complexes, defect it states are 
partially occupied. In such cases, the possibility of spin polarization may be an interesting 
issue. Some detailed discussion is given about it. We also present discussions on the elec- 
tronic structures based on the local bond picture and try to shed light on the origin of defect 
states from a different viewpoint. Recently some theoretical works have been published on 
the detailed electronic structures of MV and N doping in grapheneJ22H23l We make comments 
on their results at relevant parts in the later sections. 

Scanning tunneling microscopy (STM) is a powerful experimental technique to investigate 
the electronic features of graphene induced by individual NPD and N dopantP^ES! w ith an 
atomic resolution. In particular, the bright spots in the STM image under a small negative 
bias voltage (e.g., -0.1 or -0.2 V) would directly provide useful information for the possible 
catalytically active sites in N-graphene for ORR. Although the STM images for the clusters 
of two or more N dopants were reported, the atomic configurations were not yet solved.^ In 
addition, the STM image simulations for the configurations of an isolated graphitelike N and 
multiple pyridinelike N dopants at vacancies were also reported. 1 27 1 28 1 The actual distribution 
of different N bonding configurations in N-graphene in experiment would depend strongly on 
the energetic stabilities. In the present study, based on our results for the energetic stabilities 
of N dopants in graphene,^ we investigate the STM images of different N bonding in the 
most stable configurations including the pyridinelike N at vacancies, the pyridiniumlike N 
at MV, and the modified-graphitelike (m-graphitelike) N at a five-membered ring of 5-8-5 
DV and SW defect. 

The spectroscopy measurements of the N Is core state have indicated that N atoms are 
incorporated into graphene in the form of pyridinelike N, pyridiniumlike N, pyrrolelike N, 
and graphitelike NJ 2 * 4 * 25 * 29 -^ Note that the first two configurations of doped N would exist 
at the edges or the defect sites only. The possible N configurations at the edges of graphene 
nanoribbons and clusters have been investigated before to understand the effect of edge 
states on the N doping^ESI an( j the electron transport properties.^ In our previous worlP^ 
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as well as in our recent one we have theoretically analyzed the N Is x-ray photoelectron, 
absorption and emission spectra (XPS, XAS and XES) of N dopants along edges of graphene 
nanoribbons. In this work, we extend our analysis for the Is XPS and XAS spectra of N 
dopants near NPDs in the bulk graphene. 

The remainder of this paper is organized as follows. In Sec. [TTJ we introduce the computa- 
tional methods for electronic structure calculations and for the simulations of STM images 
and the N Is XPS and XAS spectra. Our computed electronic structures, STM images, 



and N Is XPS and XAS spectra of N-graphene are presented in Sec. |III| Finally, we draw 
conclusions in Sec. IIV1 



II. METHOD AND COMPUTATIONAL DETAILS 

In the present work, the electronic structure calculations and STM image simulations 
have been performed with the PWSCF code of the Quantum ESPRESSO suited The de- 
tailed computational setup was given in our previous paperP We only mention here that 
a 9 x 9 supercell is used for most of the calculations. The simulated constant-height STM 
images are obtained using the Tersoff-Hamann approximation.^ For a positive (negative) 
bias voltage H, the STM image provides the information for the unoccupied (occupied) 
states of the sample. A fixed sample-tip distance of d = 2 A is used throughout our STM 
image simulations. 

The calculations of the N Is XPS and XAS spectra are performed with the CP2K 
coc l e H2HSI Here a 12 x 12 supercell with a vacuum thickness of 20 A is employed in the XAS 
calculations to avoid the interference between the electrons excited to unoccupied states as 
much as possible. The atomic structures of N-graphene are taken from those optimized with 
the PWSCF codeP Two different ways are employed to estimate the N Is binding energy 
[-Eb(N Is)] for XPS spectra: (i) E^(N Is) is obtained by the Kohn-Sham (KS) energy level 
of ground state for N dopants to take into account the initial state effect and in the present 
study it is given with respect to the Fermi level; (ii) -Eb(N Is) is calculated from the total 
energy difference (ASCF) between the excited states with a full core-hole and the ground 
state and then a correction is taken into account to change the reference from vacuum to 
the Fermi level. We make some further comments on the method of ASCF calculation in 
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the present work. The core-level binding energy with respect to the Fermi level is given by 



E h (N Is) = E tot (n c -l,n v )- E tot (n c , n v ) -W, (1) 

where the first term in the right-hand side of Eq. ([T]) is the total energy of a system with 
one core-hole, the second term is the total energy of the unperturbed ground state with 
n c (n v ) denoting the number of core (valence) electrons, and the last term W denotes 
the workfunction. In our previous workj^ we evaluated W as the difference between the 
Hartree potential at the middle of the vacuum region of the supercell and the Fermi level 
of the unperturbed system. If the size of supercell is large enough, this approach does not 
have ambiguity. However, we have found that with the current supercell size, the actual 
value of the core-level binding energy depends on the size of the vacuum region. One of 
the reasons for this inconvenient aspect of Eq. ([!]) may be due to the different numbers of 
electrons in the first term and the second term, though a uniform negative background is 
added to compensate the removed core electron. To avoid this inconvenience, we used the 
following expression for W: 

W = E tot (n c , n v -l)- E tot (n c , n v ), (2) 

where the first term on the right-hand side of Eq. ^ is the total energy of a system with one 
electron removed from the Fermi level {E-p). Then, the core- level binding energy is simply 
given by 

E h (N Is) = E tot (n c - 1, n v ) - E tot (n c , n v -l). (3) 

This expression is simple and easy to be calculated. Both systems appearing on the right- 
hand side of Eq. (|3j) have the same number of electrons, i.e., N — I and we have found 
that the supercell size dependence is virtually removed. Other computational procedure 
and numerical algorithms are basically identical to those in our previous publication.^ 

III. RESULTS AND DISCUSSION 
A. N doping of defect-free graphene 

First, we examine as a reference how the atomic and electronic structures of perfect 
graphene are altered when graphitelike N is introduced. Our calculations predict the N-C 
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bond length of 1.41 A, which is slightly shorter than the C-C bond length of 1.42 A in perfect 
graphene partly due to the smaller atomic size of N. Figure [T] shows the electronic structures 
of perfect graphene and N-graphene. For perfect graphene, the original K and K' points in 
the first Brillouin zone (BZ) of primitive cell are folded into the BZ center (i.e., V point) of 
9x9 super cell, which is used in our calculations. It is known that for perfect graphene Ep 
coincides with the Dirac point, and thus both the 7r and the ir* bands near Ep are doubly 
degenerate [see Fig. fija)]. For N-graphene, Ep shifts up into the conduction band due to 
the electron donation from doped N to the ir* bands of grapheneP^I The doped graphitelike 
N also removes partly the degeneracy of the ir and tc* bands near the Dirac point. The 
impurity resonant state appears at about 0.4 eV, hybridizes with one of the bulk tt* bands, 
and forms two flat bands at about 0.18 and 0.59 eV above Ep, as shown in Figs. []Jb) and 
[jjd). The supercell size dependences of these flat bands were studied using 3n x 3n supercell 
with n changing from 2 to 5. We found that the upper flat band energy Ei and the lower 
one Ei are very well fitted (with the least-squares fitting) by 

E 2 = 2.416/n + 0.115, 
Ei = 1.264/ra + 0.119, 

where both energies are measured in eV with reference to the Dirac point. These equations 
suggest that in the limit of large n both energies converge to the same value given by about 
0.12 eV. Therefore, the two split peaks in Fig. [ltd) just above the Fermi level will be a 
single peak at 0.12 eV above the Dirac point in the dilute limit of doped N. The result is 
consistent with the recent STS measurement in which the resonant peak is located at about 
0.14 eV™>Tlie strong size dependence of E\ and E% suggests strong long range interaction 
among defectsp223 



B. Monovacancy and N doping 

1. MV in undoped graphene 

a. Band calculations For MV in undoped graphene, as three electrons are accommo- 
dated to three dangling a orbitals, the three carbon atoms next to the vacancy site would 
undergo a Jahn- Teller (JT) distortion, i.e., two of them [labeled C5 and C5' in Fig. |2|a)] 
will form a new bond.^^^ As shown in Ref. |22f the magnitude of JT distortion for MV 

7 



depends on the employed supercell size. Generally, a smaller supercell predicts weaker JT 
distortion. In the present study, the spin-polarized GGA calculations on a MV in a 9 x 9 
supercell predict that the C5-C5' bond length in the reconstructed MV is 1.96 A and that 
the CI atom stays in the graphene plane, which agree with the more recent calculations.^^ 
Judging from the results for larger supercells, the converged value of C5-C5' bond length 
may be 1 or 2% smaller than our valuePFi gure[2^c) presents the band structure of undoped 
graphene with MV. It can be seen that Ep is below the Dirac point, indicating that MV 
acts as a hole dopant. The double degeneracy of the 7r and ir* bands near the Dirac point 
is also removed by MV in the supercell calculations like the present one. 

The electronic structure obtained with the constraint of non spin polarization (for the 
structure optimized in the spin-polarized state) gives -4.0 and 2.8 eV with reference to Ep 
as the energy of the bonding state among the three dangling o orbitals and that of an anti- 
bonding state between C5 and C5', respectively. These states are far away from E-p and 
do not play any roles in transport and chemical activities. Only a dangling bond state at 
the remaining atom labeled CI has its energy rather near Ep which splits into an occupied 
majority spin state at -0.6 eV and the unoccupied counterpart at 1.8 eV with respect to 
Ep. While the exchange splitting is about 2.4 eV for the dangling a bond states, it is only 
about 0.2 eV in the 7r defect levels, which hybridize strongly with one branch of two 7r bulk 
bands in each spin state. The exchange splitting in dangling a bond states produces the 
magnetic moment of 1.0 /ib, while the defect ir states contribute to a non- integer value of 
0.46 fiB due to the hybridization with the bulk 7r bands. The net magnetic moment of 1.46 
/ie is in good agreement with the previous studies.^ 1 * 52 * 54 ! (Note, however, the comments on 



the 7r state spin polarization at the end of this section and the discussion in Sec HIE) 

The most important point for MV is that the center of gravity of the 7r defect states is 
located below the Dirac point. Therefore, the total number of electrons accommodated in 
the defect related levels is nearly five (three in a and about two in 7r) being larger by about 
one than four needed to charge compensation at MV. This is the reason why MV acts as a 
hole dopant. In the tight-binding model based on the orthonormal basis set for the ir band 
of graphene, the ir orbital level is equal to the Dirac point energy. Because of this, one may 
naively think that the defect ir states will have the center of gravity at the Dirac point. In 
the model calculations for MVp^ MV is regarded as an impurity with a strong repulsive 
potential U . In this type of treatment, the location of MV defect tt state is just below 
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the Dirac point for finite strength of Uq simply due to the hybridization repulsion between 
the impurity orbital and its surrounding orbitals, the latter being the defect ir states. This 
mechanism of locating the defect 7r level below the Dirac point may be correct or may be an 
artifact of the model. Below we propose another mechanism in which the orthogonalization 
effect will lead to such a lowering of the defect 7r level. 

In any ab initio method with localized basis set, the basis orbitals are nonorthogonal and 
we solve the following eigenvalue problem. 

HA = ESA, (5) 

where H is Hamiltonian expressed in terms of nonorthogonal basis, S is the overlap integral 
between nonorthogonal basis functions and A is the eigenvector. By performing Lowdin 
orthogonalization, we obtain 

S-^HS-^S^A = ES l ' 2 A, 

S = 1 + SS, (6) 
5-1/2 ~ i _ (i/2)8S. 

Therefore, the effective atomic level in the orthogonal tight-binding model which should be 
equal to the Dirac point energy Ed is given by 

E D = (S-^HS-^h « Hu - (1/2) ^(SSijHji + HySSp) = E a + 5E a , (7) 

where the subscripts i, j denote atomic sites, E a the atomic level given by Ha and SE a the 
correction term coming from nonorthogonality. It is important to note that the correction 
term 5E a is positive because the Hamiltonian matrix elements and overlap matrix elements 
have opposite sign. In the presence of a vacancy at site 0, the effective atomic level at the 
surrounding site n is given by 

= Ed + (i/2)(5S n0 H 0n + HoJS nQ ). (8) 

As the second term in the above equation must be negative, the effective atomic level at the 
surrounding sites of MV is lower than Ejy. In this way, we can naturally explain why the 
center of gravity of defect 7r states is below the Dirac point. 

The orbital decomposed partial densities of states (PDOS) of MV are shown in Fig. [2^e). 
The dangling a state near E F is strongly localized at CI, while the defect n states have 
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significant amplitude at CI, C3 (not shown), C5, C5', and C7. This is consistent with 
the tight-binding model analysis results that a single impurity in sublattice A induces an 
impurity state mostly localized in sublattice B and vice versa due to the existence of two 
nonequivalent Dirac pointsP^ The same behavior has also been found in the case of single 
substitutional N in defect-free graphene^ and in the results discussed below for other NPDs 
even after N doping. 

Since CI next to the vacancy site has dangling bond, it can be terminated by hydroge n 121 * 54 ! 
or substituted by nitrogen. 1 21 1 28 1 When CI is terminated by a hydrogen atom [see Fig. (2^b)], 
the atom CI moves out of the atomic plane by 0.89 A, the unoccupied minority spin dangling 
cr state of CI is pulled down below the Fermi level to accommodate the additional electron 
from the hydrogen. The doubly occupied a level with small admixture of n states is now 
merged into the bulk a band and the a spin polarization observed for non-hydrogenated 
MV disappears.^ However, the sharp defect tx states mixed with a states at the Fermi level 
can spontaneously produce spin polarization, whose magnitude is estimated to be 0.55 [1b in 
the present calculation. Nevertheless, as the energy gain from the non-spin-polarized state 
to the spin-polarized one is only about 5 meV/MV, the spin polarization will disappear at 
room temperature. (Note also the comments at the end of this section.) The band structure 
of the hydrogenated MV (H-MV) and its PDOS are shown in Figs. |2](d) and[2^f). Except the 
position of the a level and exchange splitting, the band structure of H-MV [Fig. W[d)] is very 
similar to that of MV [Fig. (2](c)] and the filling of 7r bands is not modified by hydrogenation. 
However, the peaks near Ep in the PDOS of CI have some components of a states for H-MV 
because of the out-of-plane displacement of CI. 

We make some additional comments on the spin state of MV. It is claimed in a recent 
workpZI that the spin polarization of the 7r state will vanish in the dilute limit of vacancy 
concentration. As the defect it state strongly hybridizes with the bulk 7r band, its tail must 
be rather extended. This suggests the need of careful treatment of magnetic properties of 
defect 7T states in the dilute limit. Moreover, as mentioned above in the context of H-MV, 
even if the defect tt states may be spin polarized, the magnetic stabilization energy is only 
on the order of meV. We also find that when the hydrogen atom and its nearest-neighbor 
CI atoms are confined in the atomic plane of graphene, the a states of the CI atom do 
not mix with the sharp 7r states near Ep and the spin polarization of the it states of H- 



MV disappears. Therefore, hereafter up to Sec. |IIip, we do not take into account the spin 
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polarization in the defect tt states in similar situations. The problems of spin polarization 



are summarized in Sec. Ill E for all the subtle cases treated in the present work. 

To examine the localized states near Ep, the simulated STM images under different 
bias voltages are presented in Fig. [3j They serve as a powerful test of the calculated results. 
Occupied (unoccupied) 7r states are probed by negative (positive) bias voltage. Unless we use 
a very small tip-sample distance, a states which are mostly confined within a graphene plane 
are not visible by STM. Therefore the STM images reflect the 7r state components of PDOS 
shown in Figs. (2[e) and[2^f). Figure [3]^a) shows the simulated STM images for the undoped 
MV. The pattern in the STM images under both the positive and the negative bias voltages 
exhibits a threefold symmetry with respect to the vacancy site and significantly bright 
spots at C3, C5, and C7 atoms. It is noted that the brightness in the STM images under 
the negative bias voltages is slightly stronger than that under the positive bias voltages, 
indicating p-type doping behavior for the undoped MV. The simulated STM images for the 
undoped H-MV are presented in Figs. |3](b) and^c). As the hydrogenated carbon atom CI 
is out of the atomic plane, we study two cases with the tip above and below the atomic 
plane with the same distance of 2 A in the STM simulations. In the case of the tip above 
the atomic plane in the STM simulations, the CI atom and the attached hydrogen atom are 
at the same side with the tip, and thus the CI atom shows a very strong bright spot, while 
no bright spot appears at the CI site in the case of the tip placed blow the atomic plane. 

b. Local bond picture Although it defect states are fairly extended, we demonstrate that 
the local bond picture which gives us insight to the electronic structures from a different 
viewpoint is quite useful. 

Figure [4] shows five possible ways of arranging Clar sextets^^ and double bonds in 
the graphene network with MV. As the system has a mirror symmetry with respect to the 
line AA', there is another counter part connected by the mirror symmetry for each pattern 
shown in Figs. [l|b), [l|c), |4](d) , and |4](e) . In these figures, a blue lobe denotes a singly 
occupied dangling a orbital attached to CI. The dashed circle means that the Clar sextet is 
somewhat distorted (called pseudo Clar sextet hereafter) and the triangle shows the atom 
with a dangling 7r orbital. We can see that the sites with the triangle correspond to those 
with partially filled sharp 7r DOS which are bright in the STM image with a bias voltage of 
-0.2 V. Moreover, the triangles exist at C5 and C5' in both Figs.|4]^b) and|4]^c) and the weight 
of the configuration, Figs. |i](d) and|i](e), will be small because of the lower probability for 
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the double-bond formation at the elongated C4'-C5' bond (1.49 A). These observations are 
consistent with the strongly highlighted spots in the STM image. 

The dangling bonds discussed here are the consequence of a given atomic structure. 
Therefore, once a new chemical bond with another atom or molecule is formed at a given 
dangling bond, both the electronic and the atomic structures are significantly modified. 
However, we expect that the presence of a dangling bond will contribute to the lowering of 
the activation barrier (if any) when a foreign atom or molecule approaches the atom with 
a dangling bond. In this sense, the distribution of dangling bonds predicted by the local 
bond picture gives us useful information about the chemical reactivity. The Clar sextet 
patterns of Fig. [4] imply that although a dangling bond may exist at different sites, there is 
no possibility of simultaneous existence of more than one dangling bond. This means that 
although both C5 and C5' sites are bright in the STM image for a bias voltage of -0.2 V, 
they are not simultaneously chemically active. Because of this, a diatomic molecule like O2 
will approach the MV with one oxygen atom heading for either C5 or C5'. More detailed 
analysis on the oxygen molecule adsorption at MV will be presented in a future publication. 

If a CI atom is hydrogenated, simply the blue lobe is doubly occupied to form a lone pair 
and can be removed from Fig. |4| Nothing happens in the 7r states. 

2. N doping at MV 

a. Band calculations Now, we discuss the electronic structures of N-graphene with MV. 
The results presented below correspond to the most stable configurations of substitutional 
N dopants near MV. 

The band structures of graphene doped with one and two substitutional N atoms at 
MV are presented in Figs. [5]^d) and[5|e), respectively. When one of three C atoms next 
to the vacancy site (e.g., CI) is substituted by N to form a pyridinehke N [see Fig. |5[a)], 
the unoccupied minority spin dangling a bond state is pulled down below E-p to make the 
er state nonmagnetic similarly to the case of H-MV. The C5-C5' bond length is reduced 
by 0.18 A. However unlike H-MV, N stays in the plane and the a dangling bond state is 
located in the bulk it band region to form a dispersionless sharp level rather close to Ep. 
Except these differences in the a states, the hydrogenated CI and the pyridinehke N behave 
quite similarly for 7r states and the filling of ir bands is basically the same as that of MV. 
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Interestingly, the bonding 7r bands are completely filled by two N dopants (i.e., Nci+Nc7). 
For H-MV, N doping [pyridiniumlike N, Nqi, as shown in Fig. [5^c)] fills the bonding n bands 
similarly to the two N dopants at MV. 



While the PDOSs corresponding to the above three band structures are given in Fig.Sl 
in the Supplemental Material,®' we show the simulated STM images for N-doped MV and 
its hydrogen terminated case in Figs. [61(a) and [6^b). For the N-doped MV, C5 and C7 
atoms are bright for V& = —0.2 V. The problems of 7r-state spin polarization in this case 



are discussed in Sec. HIE. On the other hand, the STM images for pyridiniumlike N at 
MV under = —0.2 eV and Vb = +0.2 eV show much lower contrast reflecting the band 
structure of Fig. |5](f). The situation is the same for the two N dopant case of Figs, |5^b) and 
Be). 



b. Local bond picture How can the situations described above be explained with the 
local bond picture? For N-doped MV, we simply remove the blue lobe in Fig. [4] like for H- 
MV and there are no changes in the distribution of Clar sextet circles and double bonds. For 
the case of two N dopants shown in Fig. [5|b), we show two distributions of Clar sextets and 
double bonds in Fig. [7| In both configurations, there are no dangling bonds being consistent 
with the band structure of Fig. ^e). However, compared with Fig. [7^ a), one pseudo-Clar 
sextet circle is missing and several double bonds are fixed in Fig. Ujb). Therefore, we 
expect that the configuration of Fig. [7]^a) may have a larger weight in the ground state. The 
situation is slightly different for pyridiniumlike N at MV [Fig.[5](c)]. The configuration shown 
in Fig. [8^a) and its counterpart obtained by mirror symmetry do not have any dangling bonds 
being consistent with the band structure of Fig. |5^f). However, C4'-C5' distance is elongated 
(1.48 A) and the probability of double bond formation at this bond will be reduced. On the 
other hand, the configurations of Figs. [8](b) and [8|c) have two dangling bonds and one less 
(pseudo) Clar sextet. These arguments suggest that pyridiniumlike N at MV is less stable 
than the case of two N dopants. In fact, the defect 7r states are located closer to Ep for 
pyridiniumlike N than for two N dopants. 
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C. Divacancy and N doping 



1. DV in undoped graphene 

The ordinary configuration of DVs in graphene which contains two pentagons and one 
octagon [i.e., 5-8-5 pattern as shown in Fig. W[ a)] can reconstruct into more complex ones 
such as the 555-777 and 5555-6-7777 patterns. 1 21 1 61 1 The stabilities of these patterns of DVs 
before and after N doping have been discussed in our previous paper.^ Although the 555- 
777 DV is the most stable one among the different patterns of DVs, N doping can enhance 
the stability of 5-8-5 DV against the 555-777 one. Therefore, we treat only 5-8-5 DV in the 
present work. 

a. Band calculations C2 and C2' atoms next to DV form a new bond with a length 
of 1.73 A [see Fig. [9^ a)]. The band structure and PDOS for graphene with 5-8-5 DV are 
presented in Figs, ^b) and |9](c). As for the a states, there are two sharp peaks in the 
PDOS of C2 at about 4 to 5 eV above Ep. These states have anti-bonding character 
between C2 and C2' of Fig. [9]^a) and their bonding counterparts are located in the energy 
range of the occupied bulk a band to become resonance states. Therefore four electrons are 
accommodated to the lower two defect resonance a states, so that the charge consistency 
is satisfied within a states. A rather subtle feature exists for tt states. It is interesting to 
find that a quite flat defect tt band exists just above Ep and that no carriers are created in 
the bulk ii bands. As CI, C2, and C4 belong to the same sublattice, the p z states of CI, 
C2, C4, and their equivalent ones are the dominant contribution to the flat defect tt band. 
In Ref. E21 however, the authors found that Ep is below the Dirac point and crosses the tt 
bands of graphene. However such results were obtained by a smaller (i.e., 7x7) supercell 
and a localized orbital basis. The use of a larger supercell and a plane-wave basis makes our 
result more reliable. 

We make two comments on this flat defect tt band. First, the flatness of the band 
does not necessarily mean the defect state is very localized. This defect state strongly 
hybridizes with one of the bulk tt bands, whose modification is a measure of the strength 
of hybridization. Therefore, like in the case of MV, the tail of the defect tt state may be 
rather extended. Second, the meaning of the presence of this flat defect tt state just above 
E F can be understood in the following way. The defect tt state just above Ep has the same 
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symmetry as that of the bonding tt molecular orbital of the dimer removed from DV as 
suggested by the STM image with V h = +0.2 V [see Fig. |9](d)]. The implication of the flat 
defect ii state being unoccupied is that the two n electrons originally belonging to the dimer 
are also removed. Therefore, the charge consistency is satisfied also within tt states without 
shifting the Fermi level with regard to the bulk tt bands leading to no carrier formation. 
This is distinctly different from the case of MV where the defect 7r level is mostly below Ep. 

Figure |9^d) shows the simulated STM images for some different bias voltages as indicated 
in each panel. As the defect level is not occupied, the STM images with negative bias voltages 
have much weaker contrast than those with positive bias voltages. 



b. Local bond picture Figure 10 shows the four ways of Clar sextets and double bonds 
around 5-8-5 DV. For the pattern in Fig. [lO^a), there is no dangling n orbital by considering 
a double bonds for CF-C2' and Cl"-C2". However, as these bonds (about 1.47 A in length) 
are elongated, the pattern may not be strongly stabilized. On the other hand, configurations 



shown in Figs. 10 'b)-(d) do not have these double bonds but have two sites with dangling tt 
states near the vacancy. Therefore we expect that these three patterns may also mix with 



the one in Fig. 10 a) and they may correspond to the presence of defect states just above 
the Fermi level. 



2. N doping at DV 

a. Band calculations Single substitutional N dopant at these DVs prefers the vertex 
site of a pentagonal ring to form a m-graphitelike N bonding configuration as shown in 
Fig. [TT|a). We consider also a pyridinelike N of Fig. [TTjfo) to analyze the site dependence. 
In particular, two and four pyridinelike N atoms at the 5-8-5 DV are favorable to aggregate.^ 
Therefore, it is interesting to see how the electronic structures of the 5-8-5 DV changes with 
N doping. The band structures for one, two, and four substitutional N atoms at the 5-8-5 
DV are presented in Figs. fTTVe) - (h). 



For a single substitutional N dopant at C4 [i.e., N C4 , Fig. 11 a)], the defect band of 5-8-5 
DV is pulled down and becomes more than half filled [see Fig. 11 'e)], producing holes in 
the bulk ii bands. Therefore, interestingly a single m-graphitelike N at 5-8-5 DV acts as an 
acceptor rather than a donor. For Nc4 at 5-8-5 DV, the C2-C2' bond length is shortened 
by 0.11 A, and thus the corresponding unoccupied a states of this C-C pair are pushed up 
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toward higher energy [see the PDOSs in Fig. S2(a) in the Supplemental MateriaP^]. For a 
single pyridinelike N at 5-8-5 DV [i.e., N C 2, Fig. [TTj^b)], it is less stable than the configuration 
of N C4 by about 0.97 eVp] A single pyridinelike N at 5-8-5 DV does not possess a lone-pair 
of a electrons, in contrast to the one at MV. Nc2 still forms a a bond with C2'. However 
the unoccupied a states of this Nc2-C2' pair shift toward lower energy (see the PDOSs in 
the Fig. S2(b) in the Supplemental Material) due mostly to a deeper potential of N. The 
two band structures [Figs. [TT](e) and[TTjT)] look similar near Ep. A single pyridinelike N at 
5-8-5 DV also acts as an acceptor. As the defect 7r states are partially occupied for both 
N C2 and N C4 , the problem of 7r-state spin polarization is discussed in Sec. |HIp ). 

For two N dopants at 5-8-5 DV, they prefer to form a dimerized pyridinelike N config- 



uration as shown in Fig. 11 [c). We first note in Fig. [lT[g) that two flat o bands exist at 
-0.4 eV and -2.3 eV, which correspond to antibonding and bonding states between dangling 
a orbitals at N C2 and N C2 '. Because both bonding and antibonding states are occupied 
for two N dopants, the a bond between them does not contribute to their binding and 
the Nc2-Nc2' distance (2.51 A) is longer than the original C2-C2' distance by 0.78 A and 
even longer than the unrelaxed C2-C2' bond of 2.46 A. From the configuration of a single 
pyridinelike N [Fig. |TTjfo)] to that of a dimerized pyridinelike N [Fig. 11 'c)], the antibonding 



a band is pulled down below Ep, so that two o electrons are added. To compensate for 
the over-screening, the defect 7r band is pushed up to remove one 7r electron and the Fermi 
level moves back to the Dirac point. A nearly flat defect band appears at 0.2 eV above Ep 
and it comes mainly from the contribution of the p z states of the dimerized pyridinelike N 
atoms and the dimerized C atoms next to DV. This can be seen from the simulated STM 
images under different bias voltages. As shown in Fig. [TT|(i), only the C2"-C2 /// pair at the 
five-membered ring looks bright by applying positive bias voltage of 0.2 V. 

The electronic structures for tetramized pyridinelike N at 5-8-5 DV [Fig. [TT^d)] are similar 
to those of dimerized pyridinelike N at 5-8-5 DV. The N-N bond becomes even longer to 



be 2.64 A. The four flat defect a states in Fig. 11 ^h) correspond to the four bonding and 
anti-bonding combinations among four dangling a orbitals. From Fig. [TTj^g) to Fig. [Ti~|(h), 
one more antibonding a state is occupied and two additional electrons introduced by two 
additional N dopants are accommodated there. In both cases of dimerized and tetramized 
pyridinelike N at 5-8-5 DV, no carriers are created by the N doping and the presence of the 
Fermi level at the Dirac point guarantees the stability of these configurations. The PDOSs 

16 



for N doped DV are given in Fig. S2(c) in the Supplemental Material. 

The simulated STM images for the tetramized pyridinelike N at 5-8-5 DV are presented in 
Figs. |TT](j) . Under the applied bias voltages of Vb = —0.2 V and Vb = +0.2 V, the simulated 
STM images show only the background because of no localized states in this energy range. 
When a more negative bias voltage (e.g., -0.6 V) is applied and a sample-tip distance of 
d = 2 A is employed, the pyridinelike N sites would have very strong brightness due to the 
unpaired a states. When the sample-tip distance increases to 2.5 A, those a states will not 
be visible in the STM image, because their spatial distribution is confined mostly within the 
graphene plane. 

b. Local bond picture The analysis of stability of N-doped DV and STM images based 
on local bond picture are given in the Supplemental Material (Fig. S3 for the simulated 
STM images of a m-graphitelike N at 5-8-5 DV; Fig. S4 and S5 for the local bond pictures 
of a single m-graphitelike N and dimerized pyridinelike N at 5-8-5 DV, respectively.)^ for 
Figs. [TTj^a) and[TT^c). We present the local bond analysis only for the tetramized pyridinlike 
N of Fig. [njd). 

Figure [12] shows two ways of arranging the Clar sextets and double bonds around 
tetramized pyridinelike N at 5-8-5 DV. The comparison of this case with the undoped DV 
clarifies the characteristic features of the tetramized pyridinelike N configuration. In the 
present case, the N-N bond corresponding to the C2-C2' bond in Fig. [9ja) is broken, which 
makes the Nc2-Cl distance shorter (1.325 A). Therefore, the configuration of Fig. 
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a is 



quite stable. Moreover, the pattern in Fig. |12|(a) has no dangling bonds. Therefore, this 



pattern is more stable than that in Fig. 12 b). 



D. Stone- Wales defect and N doping 



The SW defect is formed by in-plane 90° rotation of the Cl-Cl' bond [see Fig. 13 a)]. 
This rotation not only changes the bonding direction of CI to its nearest neighboring C2, 
but also results in the shortening of the Cl-Cl' bond by 0.10 A and the elongation of C1-C2 
bond by 0.04 A. As seen from the band structure of graphene with the SW defect [Fig. [f^d)], 
the Fermi level coincides with the Dirac point and no additional carriers are introduced into 
graphene]^ Nevertheless, the SW defect also induces a defect ir state which is located about 
0.4 eV above E F and hybridizes with one of the unoccupied bulk n bands. The defect n state 
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is localized significantly at the CI and CI' sites as seen from the simulated STM images 



for different bias voltages in Fig. 13 g) and more clearly in PDOS given in Fig. S6 in the 
Supplemental Material. Under the bias voltages of ± 0.2 and —0.5 V, the simulated STM 
images show only the background, indicating that no localized defect states appear below 
Ep + 0.2 eV. Under the applied bias voltage of +0.5 V, the bright spots in the simulated STM 
images form a pattern with rectangular shape and extend along the bond direction of the 
rotated C-C dimer. In particular, the rotated carbon atoms have strongest brightness in the 
simulated STM image. Two ways of Clar sextets and double bonds around SW defect are 



presented in Fig. 14 There are no dangling ti orbitals in both ways. This may correspond 
to no defect states below the Fermi level, as discussed above in the band calculations. 
Considering the Cl-Cl' bond is significantly shorten after rotation in SW defect, a double 
bond may exist at the Cl-Cl' dimer. Therefore the pattern in Fig |l4|a) is more probable 
and it may correspond to the ground state of SW defect. 

Upon N doping as in Fig. [i~3~|(b), the defect tt state comes down and partially occupied 



and small amount of electrons are doped into the bulk it band [see Fig. 13 e)]. Therefore, a 
single m-graphitelike N at SW defect acts as an electron donor. From the simulated STM 
images shown in Fig. [l3)4i), we can see that the CI, CI' and C2 sites are brighter under 
the bias voltage of ±0.2 V, suggesting these carbon sites may be more chemically reactive 
after N doping. The possibility of spin polarization of the defect it state in this case will be 
discussed in Sec. |III[E. The defect band is pushed completely down in close proximity to Ep 
and the possibility of spin polarization disappears after doping two N atoms [see Figs. [l~3~|(c) 
and (f)]. 

From the STM images and the PDOS given in Fig. S6 in the Supplemental Material, 
we see that the DOS at Ep is large in the order of CI', C2 and CI for a single m-graphitelike 



N at SW defect in Fig. 13 b). As there are many possible ways of arranging Clar sextets 
and double bonds in this case, it is difficult to give a consistent picture of explaining the 
behaviors seen in PDOS and STM images using local bond picture. Among many possible 
configurations, we only show three patterns where dangling bond exists at one of CI', C2 
and CI sites. For Figs. |l5[a) and (b), only one Clar sextet circle is missing, while three are 



missing in Fig. [i~5^c). This aspect is consistent with weaker brightness and smaller PDOS 



at Ep at CI site. On the other hand, between (a) and (b) of Fig. 15, the configuration of 



(b) seems to be more probable because of the short bond length (1.34 A) of Cl-Cl', being 
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inconsistent with the trend in STM image and PDOS. At the present stage, we do not have 
a proper explanation for this inconsistency. More subtle aspects may have to be taken into 
account. 

E. Magnetic properties of nitrogen doped at native point defects 

As briefly mentioned in Sec. |IIip , Ref. [221 studied the concentration dependence of the 
spin polarization of defect tt state for MV and claimed that the tt state spin polarization 
would vanish in the dilute limit. The MV case has subtle aspects. As shown in Figs. [2^ a), 
[2|c), and [2](e) , the dangling a state is definitely spin polarized to have 1.0 /ie, which is 
maintained even in the dilute limit because of the strong localization of defect a state. 
The a spin polarization produces the spin dependent exchange potential which is felt also 
by the defect tt state. This means that even if the defect tt state is not spontaneously 
spin polarized, the tt state spin polarization will be induced as long as the tt state is not 
completely filled. Looking at Fig. 3 of Ref. [221 we n °t e that as the unit cell size (defect 
concentration) increases (decreases), the Fermi level increases to the Dirac point. On the 
other hand, the energy of the defect tt state also increases. Therefore, a very careful and 
rigorous study is needed to check whether the tt state polarization will vanish or not in 
the dilute limit of MV. Note that we are discussing quantities of 0(1/ N) with N denoting 
the number of atoms in the supercell. The concentration dependence of magnetic moment 
shown in Fig. 8 of Ref. |22| clearly shows that the 7r-state spin polarization decreases with 
decreasing the MV concentration. Nevertheless, the magnetization value (including the a 
contribution) in the dilute limit in the figure may not be 1.0, although the authors fitted 
the calculated data with a function which takes 1.0 at the zero concentration limit. A more 
rigorous treatment seems to be needed to reach a definite conclusion about the limiting value 
of the 7r-state spin polarization. 

The subtlety of the defect 7r-state spin polarization at MV is caused by the existence of the 
defect u state spin polarization. In this context, the H-MV case needs careful consideration, 
because the out-of-plane displacement of H and CI in the spin-polarized state in Figs. |2](b), 
[2]^d), and [2^f) produces o-ix hybridization leading to a spin polarization at the CI atom. 
As already mentioned in Sec. |III[ B, the spin polarization totally vanishes both in a and tt 
states once H and CI are confined within the atomic plane of graphene. Moreover, as the 
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weight of a component in the defect states at the Fermi level is very small, the role of a 
state is different from that in the nonhydrogenated MV case. Only by increasing the unit 
cell size from 9 x 9 to 12 x 12, the ground state of H-MV is not spin polarized (see Table [T]). 

Other cases of possible 7r-state spin polarization are MV [Figs. |5](a) and (d)], 5-8-5 DV 
[Figs. [11(a) and (e); Figs. [DJb) and (f)] and SW defect [Figs. [l3|(b) and (e)] for all of 
which single N is doped at the defect. We first note that in all of these cases, there is no 
a spin polarization. We checked the spin-polarized n states for these cases and the results 
are summarized in Table |TJ It is found that a pyridinelike N at MV does not induce spin 
polarization. This is similar to the case of H-MV with H and CI confined in the atomic 
plane of graphene. For m-graphitelike N atoms doped at 5-8-5 DV and SW defect, the ir 
defect level is spin polarized if we use the 9x9 supercell, although the spin polarization 
energy is only of the order of meV. However, the spin polarization is significantly reduced 
or totally suppressed by increasing the supercell size to 12x12. Therefore, we can safely 
conclude that the 7r state spin polarization will vanish in the dilute limit of 5-8-5 DV and 
SW concentrations. 

F. Simulated N Is core level spectra 

Finally, we discuss N Is core level spectra for different N bonding configurations near 
NPDs. 

1. N Is binding energy and its chemical shift 

Table [H] presents the chemical shifts of N Is core level binding energy for substitutional 
N dopants around NPDs with respect to graphitelike N in defect-free graphene. The results 
obtained by the Kohn-Sham (KS) energy level of the ground state take into account the 
initial state effect. On the other hand, the ones from the total energy difference (ASCF) 
based on Eq. (|3| include the final state effect. Here our results indicate that although 
the two different estimates give about 20 eV difference in the N Is binding energy, its 
chemical shifts are determined mainly by the initial state effect, consistent with the trend 
of N dopants along edges of graphene nanoribbon d 38 * 39 ! and nitrogen-containing molecules^. 
This is strikingly different from the situation in the C Is binding energies in graphene, 
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where the final state effects govern the chemical shift We do not understand quite well 
the reason of this difference. Maybe the core-state binding energy of an impurity atom and 
that of a host atom may have qualitative difference in the final state effect. It is noted that 
the N dopants at the vertex sites of a pentagonal ring (i.e. m-graphitelike) in 5-8-5 DV and 
SW defect have quite small chemical shifts with respect to the graphitelike N. This small 
chemical shift is generally smaller than the peak width of the experimental spectra. The N 
Is core level of pyridinelike N at MV and 5-8-5 DV is about 3.8 eV shallower than the one of 
graphitelike N, which qualitatively agrees with the experimental resultsP^l The chemical 
shift of pyridiniumlike N at MV is about —1.2 eV, suggesting that this type of pyridiniumlike 
N will contribute to the XPS peak generally assigned to pyrrorlike N. We noticed that the 
chemical shift of pyridiniumlike N has significant dependence on its environment: According 
to our recent workj^ the corresponding values are —1.0 along the zigzag edge and —0.2 
along the armchair edge. 



2. N K-edge XAS spectra 

For simulating XAS spectra, we used a larger unit cell of 12 x 12 to reduce the interfer- 
ence between the excited states at neighboring cells. Even with such a large unit cell, the 
interference effect cannot be totally removed as we point out in the following arguments. 



The XAS spectra are shown for tt* and a* components separately in Fig. 16 



a. Graphitelike N The simulated N ff-edge XAS spectrum for an isolated graphitelike 



N in graphene is presented in Fig. 16 a). Two main features are found in the theoretical 
spectrum, i.e., a sharp peak at 400.8 eV which is due to the 7r* resonance and an intense 
peak at 407 eV which is contributed by the a* resonance. The positions of the ir* and a* 
peaks in our calculations agree well with the experimental results.^- 5 * 33 * 34 ! Small shoulders and 
a valley around 403 eV are the result of the interference effect mentioned above. 

b. M-graphitelike N at pentagon We examine the m-graphitelike N at two different 



types of NPDs, i.e., the 5-8-5 DV and SW defect and show their XAS spectra in Figs. 16 b) 
and [l6^c). The atomic structures for such N bonding configurations are presented in 
Figs, [ill a) and [L3](b) . The main features of two spectra look very similar, including the 
position and the profile. The sharp feature at 401 eV is attributed to N Is— > 7r* and the 
broad feature at 407 eV is associated with N Is— > a*. It is noticeable that the position 
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and profile of the two main features of m-graphitelike N are also quite similar to those of 
graphitelike N, although both peaks are less sharp and the edge of o* component is shifted 
to lower energy by about 1.0 eV for m-graphitelike N. 

c. Pyridiniumlike N at MV The simulated XAS spectrum for single pyridiniumlike N 
at MV is presented in Fig. |l6^d). This looks quite similar to the one for m-graphitelike N 



at SW defect [Fig. 16 'c)] except lower energy shift by about 1.0 eV. We noticed that the 
present XAS spectrum is quite different from those along edges.®! j n ^he latter there is a 
sharp peak of N-H a* origin between 7r* and a* peaks. 

d. Pyridinelike N at vacancies We first consider three cases, i.e., single pyridinelike N 
at MV, two and four pyridinelike N at 5-8-5 DV. Their atomic structures can be referred to 
Figs. |3~|(a), [TT](c) andfTTj^d). The simulated XAS spectra are presented in Figs. 16 'e), [l6~)T), 
and [l6|(g). For a single pyridinelike N at MV, there is a tiny peak at 396.5 eV which is 
associated with the transition to the unoccupied part of the defect ir state. A sharp peak 
appears at 397.5 eV in the simulated spectrum and it is attributed to the transition to the N 
component of the bulk ir* state. However the feature associated with a* states is quite broad. 
The XAS spectra for two and four pyridinelike N at 5-8-5 DV are almost the same, reflecting 
the similar band structures for them shown in Figs. [TTj^g) and[TT](h). The broad low energy 
7i* peak comes from the transitions to the two flat 7r* bands. However, the separation of the 
two flat bands will decrease as the defect concentration decreases similarly to the case of two 
flat bands for doped N in the defect-free graphene as discussed already. Other features are 
similar to those of single pyridinelike N at MV. The position of the feature associated with 
N Is—)- 7r* for pyridinelike N in the theoretical spectrum agrees well with the one reported 
in experiment j 34 l 64 l 



IV. CONCLUSION 

The electronic structures of N-doped graphene have been studied by performing the 
DFT calculations. We analyze the doping behavior of NPDs and substitutional N dopants 
in graphene. The origin and character of the localized states induced by NPDs are also 
carefully analyzed. Our results show that MV acts as hole dopant and the hole induced 
by a MV can be compensated by a single pyridiniumlike N dopant (but not by a single 
pyridinelike N dopant) or by two N dopants. No carriers are induced by 5-8-5 DV and 
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SW defect. The doping behavior of substitutional N is affected significantly by the defect- 
induced states. Single pyridinelike N at MV and single m-graphitelike N at 5-8-5 DV both 
exhibit hole dopants, while single m-graphitelike N at SW defect acts as electron donor. 
The dimerized and tetramized pyridinelike N at 5-8-5 DV do not induce any carriers. The 
possibility of spin polarization of n state by defects and N-doped defects has been carefully 
discussed. Although some ambiguity still remains for MV for which the defect a state has 
stable spin polarization, the spin polarization of defect n state without being accompanied 
with the a spin polarization will not exist in the dilute limit of defect concentration at least 
for the defects studied in the present work. We have also discussed the characteristic features 
in the simulated STM images and the Is related XPS and XAS spectra for pyridinelike N 
at vacancies, pyridiniumlike N at MV, and m-graphitelike N at a five-membered ring of DV 
and SW defect. These will provide useful information for identifying the different N bonding 
configurations in N-graphene. Moreover, the partial density of states at various sites around 
N-doped NPDs (presented mostly in Supplemental materials) and more directly the STM 
images for small bias voltages will also give us hints about the catalytically active carbon 
sites. These results obtained by the standard first-principles electronic structure calculations 
were explained by the local bond picture, particularly in connection with STM images for 
some cases studied here. 
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TABLE I. Total magnetic moment (M t , /xb per cell) and the spin polarization energy (Esp, in 
eV) of doped nitrogen near the NPDs in the 9x9 and 12 x 12 supercells of graphene. Here we 
present the results for the pyridinelike N at MV [Nci shown in Fig. [5ja)] and 5-8-5 DV [Nc2 shown 
in Fig. \U\b)], and m-graphitelike N at 5-8-5 DV [N C 4 shown in Fig. |IT|(a)] and SW defect [N C 3 
shown in Fig. |13[b)]. The results for H-MV are also shown for comparison. E^p is defined as the 
total energy difference between the spin-polarized state and the non-spin-polarized state. 



Configuration 
N in Figure 
Defect 



N c 



5( 



pyridinelike 

N C2 , 



11 » 



N C 4, 



m-graphitelike 

N C 3: 



11' 



13 



MV 



5-8-5 DV 



5-8-5 DV 



SW defect 



H-MV 



M t , 9 x 9 
E S p, 9x9 
M t , 12 x 12 
Esp, 12 x 12 



0.00 
0.000 
0.00 
0.000 



0.21 
-0.001 
0.00 
0.000 



0.50 
-0.002 

0.00 
0.000 



0.53 
-0.003 

0.08 
-0.001 



0.55 
-0.005 
0.00 
0.000 



TABLE II. N Is core level binding energy [£7b(N Is), in eV] with respect to the Fermi level for 
graphitelike N in defect-free graphene (Nc), m-graphitelike N at 5-8-5 DV and SW defect [Nc4 
and Nc3 shown in Figs. [TT|a) and [l3^b), respectively], pyridiniumlike N at MV [Nci-H shown 
in Fig. [5^c)], and pyridinelike N at MV and 5-8-5 DV [Nci and Nc2 shown in Figs. [5^a), |Ll|c) 
and|ll[d), respectively] estimated by the initial Kohn-Sham energy level (KS) and the total energy 
difference (ASCF). The values in parentheses are the chemical shifts with respect to the graphitelike 

N. 

£b(N Is) 



Configuration graphitelike m-graphitelike 
N in Figure 

Defect defect-free 5-8-5 DV SW defect 



pyridiniumlike 



pyridinelike 



N c N C4 , 



11 



a) N C3 , 



13 



» Nci-H, 
MV 



c) Nci 



» N C2 , 



11 



c) N C2 , 



ii ;d) 



MV 



5-8-5 DV 5-8-5 DV 



KS 


379.78 


378.98 


379.17 


378.07 


375.73 


375.82 


375.83 




(0.00) 


(-0.80) 


(-0.61) 


(-1.71) 


(-4.05) 


(-3.96) 


(-3.95) 


ASCF 


400.16 


399.95 


399.93 


398.93 


396.40 


396.39 


396.41 




(0.00) 


(-0.21) 


(-0.23) 


(-1.23) 


(-3.76) 


(-3.77) 


(-3.75) 
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FIG. 1. (Color online) The band structures of (a) perfect graphene and (b) N-graphene with single 
substitutional N atom in a 9 x 9 supercell. Partial density of states (PDOS) for the 2p orbitals: (c) 
C atom in perfect graphene and (d) an isolated substitutional N atom in N-graphene. The zero of 
energy is set at E-p. Panel (b) shows the fat bands (green solid points) derived from N 2p z state. 
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FIG. 2. (Color online) The atomic structures of (a) MV and (b) hydrogenated MV [one of three 
C atoms next to MV terminated by H, denoted as H-MV]. The band structures of graphene with 
(c) MV and (d) H-MV. The partial density of states (PDOS) for the 2p orbitals of C atoms around 
defect region of (e) MV and (f) H-MV. 
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V b = -0.5 V V h = +0.5 V 

(a) 




(b) (c) 

FIG. 3. (Color online) Simulated STM images for (a) MV and (b), (c) hydrogenated MV under 
different bias voltages indicated in each panel and with a sample-tip distance of d = 2 A. In (b) 
[(c)], the hydrogenated CI and STM tip are on the same [opposite] side of graphene plane. 
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FIG. 4. (Color online) Five possible ways of arranging Clar sextets and double bonds around MV. 
The solid line AA' shows the mirror symmetry axis. The dashed circle stands for a pseudo-Clar 
sextet. The lobe and triangle symbols stand for dangling a and 7r orbitals, respectively. 
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FIG. 5. (Color online) The atomic structures of (a) a pyridinelike N (i.e., Nci) at MV, (b) a 
pyridinelike N (i.e., Nci) plus a m-graphitelike N (i.e., Nc7) at MV, and (c) a pyridiniumlike N at 
H-MV. The band structures of N-graphene with MV: (d) N C i at MV, (e) Nci + N C 7 at MV, and 
(f) Nci at H-MV. Panels (d), (e), and (f) show the fat bands (green solid points) derived from Nci 
2p z state. 




(a) (b) 

FIG. 6. (Color online) Simulated STM images for (a) pyridinelike N and (b) pyridiniumlike N at 
MV under different bias voltages indicated in each panel and with a sample-tip distance of d = 2 

A. 
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(a) (b) 

FIG. 7. (Color online) Two possible ways of arranging Clar sextets and double bonds around 
two N substitutions (a pyridinelike N plus a m-graphitlike N) at MV. The solid line AA' shows the 
mirror symmetry axis. The dashed circle and triangle symbol stand for a pseudo-Clar sextet and 
a dangling ir orbital, respectively. 




FIG. 8. (Color online) Three possible ways of arranging Clar sextets and double bonds around a 
pyridiniumlike N at MV. The solid line A A' shows the mirror symmetry axis. The dashed circle 
and triangle symbol stand for a pseudo-Clar sextet and a dangling n orbital, respectively. 
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\ 

I ,, = -(1.5 V I |, = +0.5 V \ 

(d) 

FIG. 9. (Color online) (a) The atomic structure of 5-8-5 DV. (b) The band structure of graphene 
with 5-8-5 DV. The zero of energy is set at E-p. (c) The partial density of states (PDOS) for the 
2p orbitals of C atoms, (d) Simulated STM images under different bias voltages indicated in each 
panel and with a sample-tip distance of d = 2 A. 




(a) (b) (c) (d) 

FIG. 10. (Color online) Four possible ways of arranging Clar sextets and double bonds around 
5-8-5 DV. The solid line AA' shows the mirror symmetry axis. The dashed circle and triangle 
symbol stand for a pseudo-Clar sextet and a dangling ir orbital, respectively. 
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FIG. 11. (Color online) The atomic structures of 5-8-5 DV after N doping: (a) a m-graphitelike 
N (i.e., Nc4), (b) single pyridinelike N (i.e., Nc2), (c) dimerized pyridinelike N (i.e., Nc2 + Nc2')> 
(d) tetramized pyridinelike N. The band structures of graphene with 5-8-5 DV: (e) Nc4, (f) Nc2 5 
(g) Nc2 + Nc2', and (h) tetramized pyridinelike N. The zero of energy is set at Ep. Panels (e)-(h) 
show the fat-bands derived from N 2p z and 2p x + 2p y (green solid and red open points, respectively) 
states. Simulated STM images under different bias voltages indicated in each panel and with a 
sample-tip distance of d = 2 A: (i) dimerized pyridinelike N and (j) tetramized pyridinelike N. 
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(a) * (b) 

FIG. 12. (Color online) Two possible ways of arranging Clar sextets and double bonds around 
tetramized pyridinelike N at 5-8-5 DV. The solid line AA' shows the mirror symmetry axis. The 
dashed circle and triangle symbol stand for a pseudo-Clar sextet and a dangling tt orbital, respec- 
tively. 
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FIG. 13. (Color online) The atomic structures of (a) SW defect, (b) single m-graphitelike N atom 
(i.e., Ncs) at SW defect, and (c) two m-graphitelike N atoms (i.e., Nc3+Nc3') at SW defect. The 
band structures of graphene with SW defect in different cases: (d) before N doping, (e) Nc3, and 
(f) Nc3 + Nc3'- The zero of energy is set at Ep. Panels (e) and (f) show the fat-bands (green 
solid points) derived from the N 2p z states. Simulated STM images under different bias voltages 
indicated in each panel and with a sample-tip distance of d = 2 A: (g) SW defect and (h) Nc3 at 
SW defect. 
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(a) (b) 

FIG. 14. (Color online) Two possible ways of arranging Clar sextets and double bonds around 
SW defect. The solid line AA' shows the mirror symmetry axis. 
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FIG. 16. (Color online) The computed N Is XAS spectra of substitutional N atoms in different 
bonding configurations: (a) a graphitelike N in prefect graphene; (b) a m-graphitelike N at 5-8-5 
DV [Nc4 shown in Fig. [TT^a)]; (c) a m-graphitelike N at SW defect [Nc3 shown in Fig. |l3|fo)]; (d) 
a pyridiniumlike N at MV [Nqi shown in Fig. [5j^c)]; (e) a single pyridinelike N at MV [Nci shown 
in Fig. [5^a)]; (f) dimerized pyridinelike N at 5-8-5 DV [Nc2 +Nc2' shown in Fig. |ll[c)]; and (g) 
tetramized pyridinelike N at 5-8-5 DV [shown in Fig. |ll[d)]. 
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